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We present a detailed study of soliton compression of ultra-short pulses based on phase-mismatched 
second-harmonic generation (i.e., the cascaded quadratic nonlinearity) in bulk quadratic nonlinear media. 
The single-cycle propagation equations in the temporal domain including higher-order nonlinear terms are 
presented. The balance between the quadratic (SHG) and the cubic (Kerr) nonlinearity plays a crucial role: we 
define an effective soliton number - related to the difference between the SHG and the Kerr soliton numbers 
- and show that it has to be larger than unity for successful pulse compression to take place. This requires 
that the phase mismatch be below a critical level, which is high in a material where the quadratic nonlinearity 
dominates over the cubic Kerr nonlinearity. Through extensive numerical simulations we find dimensionless 
scaling laws, expressed through the effective soliton number, which control the behaviour of the compressed 
pulses. These laws hold in the stationary regime, in which group-velocity mismatch effects are small, and 
they are similar to the ones observed for fiber soliton compressors. The numerical simulations indicate that 
clean compressed pulses below two optical cycles can be achieved in a /3-barium borate crystal at appropriate 
wavelengths, even for picosecond input pulses. © 2008 Optical Society of America 
OCIS codes: 320.5520, 320.7110, 190.5530, 190.2620, 190.4400 



1. Introduction 

Second-harmonic generation (SHG) in the limit of large 
phase mismatch gives rise to the so-called cascaded 
-^.( 2 ) ; T^(2)_ nonmlcar ^.y ; m w hich weak conversion to the 

second harmonic (SH) occurs, while a Kerr-like non- 
linear phase shift is induced on the fundamental wave 
(FW) [1,2]. An elegant theoretical view of this process is 
that in the cascading limit the system can be described 
by a nonlinear Schrodinger equation (NLSE) for the FW 
alone [3]. The induced FW nonlinear phase shift then 
comes from an effective self-phase modulation arising 
from the phase mismatch. The phase shift can be quite 
large and negative, since the phase mismatch determines 
the sign and magnitude of the effective cubic nonlinear- 
ity. Such a self-defocusing nonlinearity can be used to 
compress a pulse when combined with normal disper- 
sion [4-9], and problems normally encountered due to 
self-focusing in cubic media are avoided. Thus, having 
no power limit, a self-defocusing compressor can cre- 
ate high-energy near single-cycle fs pulses in bulk me- 
dia [4,8]. The exploited compression mechanism is the 
periodic behaviour of higher-order temporal solitons of 
the NLSE. They oscillate in temporal duration upon 
propagation, and the optimal compressor length is when 
the pulse width is the narrowest. 

As the compressor scheme exploits an effective self- 
defocusing cubic term from cascaded quadratic effects, 
the compressor will naturally be affected by the self- 
focusing cubic nonlinearity inherent to any transpar- 
ent material. This detrimental cubic nonlinearity must 
be counterbalanced and then exceeded to achieve com- 
pression [5,8]. In this work we systematically demon- 



strate theoretically and numerically how this can be ex- 
pressed conveniently using the soliton number formalism 
known from the NLSE in fiber optics [10,11]. Firstly, 
we show that the SHG soliton number [8] Ague must 
outbalance the Kerr soliton number Ak G it before com- 
pression takes place. More precisely, the effective soli- 
ton number N e g = \J -/V| HG — A£ orr must be larger than 
unity. This can be achieved by adjusting the phase mis- 
match, but only if the quadratic material nonlinearity 
is sufficiently stronger than the cubic Kerr nonlinearity. 
Secondly, when compression is successful, the compres- 
sion factor, pulse quality, and optimal compressor length 
follow certain scaling laws. Such empirical scaling laws 
expressed through Axerr were previously given for the 
NLSE [12-14]. Here we present detailed numerical simu- 
lations of pulse propagation in a commonly-used nonlin- 
ear material, /3-barium borate (BBO). We find general 
dimensionless scaling laws quite similar to the ones from 
the NLSE, except that one must express them through 
N c ff. The simulations span a very wide input-pulse pa- 
rameter space and the scaling laws therefore provide an 
important tool for experimental situations. 

A major obstacle to the cascaded-quadratic soliton 
compressor is the group-velocity mismatch (GVM) be- 
tween the FW and SH pulses. GVM induces a Raman- 
like perturbation that distorts the compression, result- 
ing in asymmetric pulses and pulse splitting [4,8,15]. 
Two distinct regimes exist: in the stationary regime the 
GVM-induced Raman-like perturbations are weak, re- 
sulting in clean compressed pulses. In the nonstationary 
regime, the GVM-induced perturbations are strong, and 
poor pulse compression is observed. Moreover, the com- 



pression behavior deviates from the NLSE scaling laws. 
Here we focus on the stationary regime, where cleanly 
compressed few-cycle pulses may be obtained, and sum- 
marize the guidelines for avoiding nonstationary behav- 
ior. Subsequent publications will focus on compression in 
the nonstationary regime in greater detail. In particular, 
we have recently devised a nonlocal theory that quan- 
titatively defines the boundary between the stationary 
and nonstationary regimes [16]. 

Finally, when describing (ultra)-short pulse propaga- 
tion in nonlinear media it is generally important to de- 
scribe both temporal and spatial effects as they will in- 
herently be interlinked. Examples include space-time fo- 
cusing in bulk media [17], and strong modification of the 
dispersion and nonlinear mode overlap in waveguide ge- 
ometries. However, the simplified ID temporal descrip- 
tion presented here is often an adequate starting point 
for understanding the general temporal behaviour of the 
system during compression. In this work we consider in- 
teraction in a bulk material (no wave guiding) and we 
neglect diffraction (i.e., the beam must not be focused 
too tightly). In subsequent publications we will address 
the behaviour in a wave-guiding geometry such as a fiber, 
as well as in bulk geometry including the spatial effects. 



scalar fields and either a Type-0 or Type-I SHG geome- 
try, in which the two FW photons have identical polar- 
ization [19]. In a bulk medium, the generalized SEWA 
propagation equations for the coupled electric fields are 
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where Ej = Ej(z,t). For a Type-0 geometry the cross- 
phase modulation (XPM) coefficient is B = 2, while for a 
Type-I geometry B = 2/3 because the FW and SH have 
orthogonal polarization [11]. We assume that the FW 
and SH spectra do not significantly overlap, an assump- 
tion generally valid down to single-cycle pulse durations. 
The linear propagation operators are 
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2. Generalized propagation equations for ultra- 
short pulses 

Below we present the generalized bulk SHG propaga- 
tion equations, including both quadratic and cubic non- 
linear effects. We aim to investigate compressed pulses 
that may be near-single optical cycle in duration, so 
we must use the propagation equations derived with- 
out approximations that directly impart constraints on 
the pulse duration and bandwidth. Such SHG equations 
have been derived elsewhere [18] using the slowly evolv- 
ing wave approximation (SEWA) [17]. Here we recast 
these equations in the framework of cascaded quadratic 
soliton compression as to define all critical experimen- 
tal parameters. We also include generalized cubic non- 
linearities from interaction of two fields with different 
frequencies (the FW and SH), see App. lAl 

The investigation is a temporal study only. We ne- 
glect transverse spatial effects, assuming that diffrac- 
tion plays a negligible role on the length scale of the 
compressor, i.e., that each point on the beam may be 
considered a plane wave. Previously-reported quadratic 
soliton compressors have lengths on the order of a few 
centimeters [4-8]. For a Gaussian beam at near-infrared 
wavelengths, this implies that our assumptions are valid 
for beam waists larger than approximately 100 fim. Sec- 
tion [F] will discuss considerations of compression with 
transverse-spatially varying profiles. Finally note that in 
our notation, a primed variable is always dimensionless. 

A. The dimensional form 

In SHG two FW photons of frequency u>\ combine to give 
a photon of frequency uj 2 = 2cji (the SH). We consider 



where Dj are the dispersion operators 
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Here kj = rijUij/c, where rij is the refractive index. The 
fields are in the frame of reference traveling with the FW 
group velocity v g> \ by the transformation r = t — z/v gt \. 
This is the origin of the GVM-term d\ 2 = l/v g .i^^/v g . 2 , 



where 
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dispersion, and Afc = k 2 — 2k\ is the phase mismatch. 
Note the unfamiliar term in Eq. (|2b[) . which is an ef- 
fective SH dispersion operator X>2,eff = V 2 + S^ 1 ^-^ 
[see Eq. (|2l)|) ]. which deviates from T> 2 due to GVM and 
self-steepening in the SEWA model [18], see App. IBl 



The quadratic nonlinear coefficient is 
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X^wi/2cnj = d e ffUJi/cnj where x is the value of the 
quadratic nonlinear tensor along the polarization direc- 
tion of the interacting waves, and as is typical we take 
d e e = x*- 2 ' /2. Dispersion of x^ is assumed negligible. 



The cubic nonlinear coefficient is 
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App. El We neglect two-photon absorption, implying 
Im(x'- 3 ' ) ) = 0, which holds as long as the SH frequency 
spectrum lies below the two-photon absorption edge of 
the medium. Note, the usual notation for riKerrj is n 2, 
but we reserve the subscript 2 for the SH. 

Self-steepening effects are modelled by the operator 
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jjt^:, It approaches unity when bandwidths 
Aujj are small compared to the carrier frequencies uij, 
i.e., for pulses longer than roughly 10 optical cycles. 



2 



Finally, TZj(t) describes the vibrational Raman re- 
sponse of the cubic nonlinearity, see App. [5] We hence- 
forth set fn = 0, because Raman scattering plays a neg- 
ligible role under the relevant experimental conditions. 

We now scale the fields by defining Aj = Ej ^eorijc/2 
so | A, | 2 gives the beam intensity Ij. Eqs. (|TJ) become 
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Kg HG = ^2u}'{d' 2 s /n{n 2 e c 3 , and j[ = Win^ err l /c. Here 
the intensity-version of the Kerr nonlinear refractive in- 
dex with unit [m 2 /W] is nj^ eII l = 2n.Kerr,i/^i£oc = 
3Re(x (3) )/4n 2 e c Note the connection to the fiber NLSE 
(where Aj is scaled to the power) ; 7/ is equivalent to the 
parameter [fO] 71 = wing err l /cA e s, where A c s is the ef- 
fective mode area of the FW. Finally, n = ni/n 2 which 
is typically close to unity, except when Ak is large. 

B. The dimensionless form 

We now rescale space and time as z' = z/Lu,i, and 
t' = r/Ti^n, where Lu,i = T?i n /\k^\ and Tu n are the 
FW dispersion length and input pulse duration. When 
compressing a FW pulse, no SH is launched, so both 
fields are scaled to the peak FW input intensity Ii t - m = 
1^(0, 0)| 2 , i.e., Uj = A 3 /y/h~. Thus, U x = £i/£i jin , 
where £ x - ln = |£'i(0, 0)|, and Eqs. become 
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where AgHG and Axem the dimensionless quadratic and 
cubic soliton numbers, respectively, are defined as 
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^Kerr = ^D,l7l^l,in = Lu,l^in^ eIll Ii; m /c, (6b) 

and the dimensionless phase mismatch is Ak' = AfcLo,i- 
The parameter Ague is deliberately chosen for a soliton 
compression scenario, as will become clear later. 

Equations {2j on dimensionless form arc C[ = 
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T>' 2 eS is given by Eq. (|27|) . The dimensionless dis- 
persion operators and corresponding coefficients are 
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and 8f l) ee LD 1 i*j m) (I^ 2 |*i 2) l m 0~ 1 Finally, using 
s' ee (u^T^in) -1 we have the steepening operators with 
dimensionless time S[ = l + is'-^p, and S 2 ee 1 -My -§p. 



C. The soliton numbers 

We now turn our attention to the soliton numbers as 
defined in Sec. EJ We will argue - and later show with 
numerical simulations - that the effective soliton number 
iVeff = y/ A| HG — A^ crr can be used to describe the sys- 
tem in general. Indeed, N c g characterizes the outcome 
for any given input state of the system. 

First, consider the NLSE in absence of quadratic 
nonlinearities. We include for instruction the Raman 
terms and ignore the steepening terms [Eq. (|25ap with 
and Sj = 1]. In dimensionless form it reads: 
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The nonlinear coefficients due to the intensity scaling are 1 q , 



.dU x sgn(/4 2) ) d 2 Ui 



Q T ,2 

|C/i| 2 f/i 



d\Ui 



0, (7) 



which holds when the pulse duration is much slower than 
the Raman response, and t' r = tr/Ti^, where t^j is a 
characteristic Raman vibrational response time [10], see 
App. [XJ The Kerr soliton number (|6b[) we express as 

Afcerr = x/^D.l/iKcrr, where LlCerr = (7l^l,in) _1 IS the 

characteristic Kerr length. Equation governs soliton 
compression in cubic nonlinear media, and the compres- 
sion characteristics are well understood to be critically 
dependent on iVxerr [10-14]. We stress that soliton com- 
pression using a pure Kerr self-focusing nonlinearity re- 

( 2) 

quires anomalous dispersion sgn(fcj ) < 0. 

Likewise, the quadratic soliton number Ague is a crit- 
ical parameter in cascaded quadratic soliton compression 
that as we shall now see has a role analogous to the cubic 
soliton number of the NLSE. Using perturbative meth- 
ods one can reduce the dimensionless SHG equations © 
- in absence of steepening and cubic nonlinear terms - 
to a single equation for the FW analogous to Eq. ([7]) 
when the phase mismatch Ak is suitably large [3, 15]. 
The approximate equation for the FW field can for soli- 
ton compression purposes be cast as [8] 
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where s\ 2 = sgn(di2). In analogy to the cubic nonlin- 
earity, a characteristic length for the cascaded quadratic 
nonlinearity Xshg = c 2 nin 2 \Ak\/uj 2 d 2 s £ 2 in , such that 
AguG = y/Ln,i/LsnG- The similarity to Eq. ([7|) is clear, 
and explains the choice of Ague in Eq. (|6ap . Similar to 
Axerr for the NLSE, Ague is namely a critical parameter 
in cascaded quadratic soliton compression. Several other 
features should be noted. First, the sign of the phase 
mismatch controls the sign of the induced cubic non- 
linearity For the purpose of soliton compression with 
normal dispersion sgn(fcj ) > 0, we require Afc > in 
order to have self-defocusing, or negative nonlinearity. 
Second, the Raman-like term in Eq. ([8]) is important in 
the compression dynamics, as will be explained in Sec.lDl 
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Equations {7J and ([8]) consider the cases where ei- 
ther the cubic or the quadratic nonlinearity is signifi- 
cant. However, we must consider both orders of nonlin- 
earity, so using the same perturbation methods, we re- 
duce Eqs. ([5]) to an approximate equation for the FW. 
To first order, and neglecting steepening effects, we find 
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In Eq. ([H]) self-phase modulation (SPM) effects from 
both quadratic and cubic nonlinear effects are present. 
For Afc > 0, x^ 2 -* : x' 2 ' self-defocusing phase shifts 
starts to cancel the x^ 3 ) self-focusing phase shifts. Thus, 
for soliton compression with normal dispersion and self- 
defocusing phase shifts, the effective soliton number 
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must be larger than unity, N e g > 1. This is analogous 
to the cubic soliton compressor, for which compression 
will only occur for A^Ken- > 1. In Sec. [3] these conclu- 
sions are supported numerically through simulations of 
Eqs. ([5]). By this we demonstrate that N e g governs all 
the behaviour of the compressed solitons, at least in the 
stationary regime. This regime is now discussed. 

D. The stationary regime 

An important observation is that Eq. §§§ has a GVM- 
induccd Raman-like term similar to the one induced by 
the delayed cubic response in the NLSE ([7]). Its charac- 
teristic dimensionless temporal response is [8, 15] 



'R.SHG 



= 2|dia|/|A*|Ti, 



(11) 



and t t? : shg = TR,SHG/^i,in. Typically, tr iS hg lies be- 
tween 1-5 fs in BBO. This term imposes a Raman-like 
red-shift of the FW, which causes the compressed pulse 
to be asymmetric and limits the compression achievable. 

To understand this better, observe the effective non- 
linear FW phase shift built up during propagation 
<I>nl(z) = -sgn(Afc) z/Lshg = -zuj 2 dl s £ 2 M /c 2 n 1 n 2 Ak 
[1], which holds when y/\Ak\LsnG 3> 1 (weak conver- 
sion to the SH). For efficient compression </>nl should at 
least be on the order of — it [4,5], and evidently the pulse 
needs to propagate a shorter distance to achieve this for 
Afc low (as long as the cascading limit is upheld). How- 
ever, Liu et al. [4] pointed out that GVM sets a lower 
limit to Afc: In phase-mismatched SHG, the phase be- 
tween the FW and SH changes several times, correspond- 
ing to a conversion/back-conversion cycle of energy be- 
tween FW and SH fields. The distance over which the 



relative phase changes sign once is exactly character- 
ized by the coherence length L co h = 7r/|Afc|. Thus, it 
is important that the temporal walk-off (GVM) length 
Lqvm = Ti,in/|di2| is (much) larger that the coherence 
length. Liu et al. [4] found that it is sufficient to demand 
that L C oh < 4Lqvm- Thereby the so-called stationary 
regime was defined, thus requiring Afc > 47r|tii2|/Ti j i n . 
Therefore it is advantageous to have low GVM because 
it gives access to large nonlinear phase shifts that oc- 
cur when Afc is not too large. In addition, low GVM 
will result in a smaller r' R SHG -tcrm in Eq. ([9]) so the 
Raman-like perturbations become smaller. This was a 
main point of Ref. [8], where it was shown that for a 
given amount of GVM the compression increases as the 
soliton number is increased until a certain critical value 
where the Raman-like perturbations start to dominate. If 
GVM is reduced the soliton number can be increased fur- 
ther before this happens, and stronger compression can 
be achieved. We recently devised a more accurate the- 
ory for the stationary and nonstationary regimes [16], 
showing that the cascaded nonlinearity gives rise to a 
nonlocal response, and that in presence of GVM the non- 
local response function is asymmetric (Raman-like). In 
the stationary regime the nonlocal response function is 
localized and Eq. ((9|) is recovered in the weakly nonlo- 
cal regime (where the response function is very narrow 
compared to the propagating pulse). In the nonstation- 
ary regime the nonlocal response function is oscillatory 
and docs not decay. This has severe consequences to the 
built up negative nonlinear phase shift and results in very 
poor compression. Finally, a more accurate requirement 
for being in the stationary regime was found as 



Afc > Afc sr = d\ 2 /2k 



(12) 



which notably is independent of Ti.in- In what follows, 
we will use this value. We should note that satisfaction 
of inequality (fT2]l docs not eliminate the GVM-induccd 
perturbation: Raman-like effects are present when GVM 
is nonzero both in the stationary regime Afc > Afc sr , 
where Eq. ([9]) holds, and in the nonstationary regime 
Afc < Afc sr . However, in the latter case, where the non- 
local response function is oscillatory and unbound, the 
Raman-like perturbations are stronger. 

E. The compression window 

In the previous sections we have imposed two require- 
ments for clean soliton compression. First, in order to 
observe solitons, the effective soliton number must be 
larger than unity N e g > N e g jC = 1. This gives from 
Eq. (p~0|) a critical value of the phase mismatch 
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Second, Afc > Afc sr to be in the stationary regime, see 
Eq. (fT2|) . which is strictly related to the size of GVM and 
SH group-velocity dispersion, and is independent of the 
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(a) T ] . =200 fs FWHM (b) I, _ =200 GW/cm 2 




I lm [GW/cm 2 ] Tj to FWHM [fs] 

Fig. 1. The compression window for BBO with Ai = 
1064 nm using (a) T 1M = 200 fs, and (b) 7i, in = 
200 GW/cm 2 . Also shown is Ak c found in numerics 
(see Sec. IA"|) . and the dashed line is an empirical scal- 
ing law (|15p based on these simulations. 

input pulse duration. This defines a compression window, 
inside which we can expect clean compression: 

Afc sr < Afc < Afc c . (14) 

Staying inside this window is therefore a question of 
choosing the right phase mismatch. Keeping the window 
open (Afc sr < Afc c ) is also a matter of having the right 
input pulse intensity and duration, see Fig. [TJ Moreover, 
there are cases where the material nonlinearity balance 
and/or GVM effects do not allow for Eq. (fl"4"|) to be sat- 
isfied. For instance, although Ak c increases with the in- 
put intensity, see Fig. [TJa), its maximum value, occur- 
ring when N^ IT -> 0, is Afc c , max = ^id 2 ff / 'c7iin 2 n K orr,i- 
Thus, the balance between quadratic and cubic nonlin- 
caritics is a fundamental limit on the peak value of Ak c . 
Moreover, when GVM is large Afc sr becomes high, so 
very large quadratic nonlincarities arc needed to open 
the window. Note that compression can be achieved out- 
side the compression window when Afc < Afc sr . However, 
due to the strong Raman-like effects in the nonstation- 
ary regime A c ff must be kept small in order to get a 
clean symmetric pulse, allowing only for moderate com- 
pression. Conversely, even inside the compression win- 
dow distorted compressed pulses may occur [16]. 

F . Transverse spatial effects 

Let us end with a brief discussion of the implications of 
including effects of modulation in the transverse spatial 
dimensions (x and y). The equations will then include 
both diffraction terms and space-time focusing terms 
(due to steepening effects of diffraction [17]). In its sim- 
plest effect, diffraction will act to alter the pulse intensity 
by altering the mode area, which will lead to changes in 
the soliton order during propagation. Therefore, when 
zq is on the order of a characteristic diffraction length 
or longer, diffraction will have significant effects on the 
compression results. Additionally, the space-time focus- 
ing terms, which couple temporal and transverse spatial 



coordinates and become significant at few-cycle pulse 
durations, will be expected to produce nontrivial spa- 
tiotemporal distortions. 

Moreover, in the plane-wave limit, i.e., where diffrac- 
tion effects arc small, while the temporal dynamics can 
be well described by the model presented in this paper, 
inhomogencity in the spatial intensity profile will result 
in inhomogencous compression across the beam. For ex- 
ample, for a Gaussian profile only the central part of the 
beam may have the appropriate intensity to be efficiently 
compressed. In order to achieve a more homogeneous 
compression, a super-Gaussian profile can be used [20]. 

Another case is where the light is guided, such as in 
an optical fiber. Here, the nonlinear terms in the propa- 
gation equations must include the effective spatial mode 
overlap areas, and thus modifying the soliton numbers 
as well. Another important difference in the guided case 
is that the dispersion is affected. Moreover, it is possible 
to achieve very high intensities with low energy pulses 
because the light can be confined in a very small area. 
We will in another publication go into details with the 
specific case of a photonic crystal fiber, and try to in- 
vestigate the implications of wave-guiding effects on the 
cascaded quadratic soliton compressor. 

3. Numerical results 

We now proceed to the numerical results, where we first 
study the transition to compression, and then derive the 
empirical scaling laws for the compressed pulses. The 
simulations were done using Eqs. ((5|), which include sclf- 
steepening, the SEWA corrected SH dispersion term and 
higher-order nonlinear mixing terms. The coupled equa- 
tions were solved using a split-step Fourier technique 
with a 2nd order Runge-Kutta algorithm to evaluate 
the nonlinear terms in the time domain. The dispersion 
terms were evaluated in the frequency domain, and since 
kj{u>) is known analytically (from the Scllmcicr equa- 
tions of BBO [19]) we can actually evaluate the disper- 
sion operator Dj exactly [21] instead of using the expan- 
sions ((3]). For the effective SH operator (|27|) we used 30 
terms in the expansion of S^.The steepening terms were 
applied by using the convolution theorem and thus let- 
ting the steepening operator act on the nonlinear terms 
in the frequency domain. The number of discretization 
points N z in the z-direction was chosen so that around 
15 steps were taken within a single coherence length 
L C oh = 7r/|Afc|. The number of temporal points were 
2 12 -2 15 , primarily dictated by the fact that due to GVM 
the SH will have a trailing pulse, which must stay inside 
the time window of the simulations. 

We decided to use BBO as the nonlinear medium for 
reasons described in App. O We chose the Yb/Nd:YAG 
wavelength Ai = 1064 nm. At this wavelength BBO has 
a medium level of GVM (d\2 — —100 fs/mm). If instead 
we consider Ai = 800 nm - another wavelength of con- 
siderable interest - GVM is roughly twice as big and the 
compression window is much smaller. From this point of 
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Fig. 2. Locating the critical transition point to pulse 
compression. The critical phase mismatch Afc c vs. A^Kerr 
in a semi-log plot for simulations having 7i,i n = 1 — 
600 GW/cm 2 and T 1M = 80 - 600 fs FWHM. 



view Ai 



1064 nm is better suited for this investigation. 



A. Transition to compression 

First, we characterize the system at the transition to 
compression, in order to understand whether the require- 
ment iVeff > iVeff, c = 1 holds. Quantitatively, we first 
found the critical phase mismatch Afc c for a given Ii t i n 
and Ti.i n (thus, A?Kerr is fixed, while Ashg is varied as Afc 
is scanned). This was done through numerics of L = 50 
mm propagation in a BBO crystal. First, we located a 
Afc just before the transition (slight decrease in pulse in- 
tensity due to pulse broadening) and a Afc just after the 
transition (slight increase in pulse intensity due to pulse 
narrowing). Afc c was then found by interpolation. We 
chose sech-shaped input pulses Ui^n = sech(T') because 
they are solitonic solutions to the NLSE (without higher 
order dispersion and nonlinearities), and the pulse in- 
tensities during propagation were correspondingly fitted 
to a sech 2 -shaped pulse. Based on Ref. [22] we expect 
the results presented here to remain largely identical for 
sech- and Gaussian-shaped input pulses. 

In Fig. [21(a) we plot Afc c vs. A?Kerr for many different 
values of Ti^n and Ii. ln . As Sec. [C] suggested, the transi- 
tion to compression can expressed as N e f{ > N e s tC = 1, 
and the corresponding theoretical line calculated from 
Eq. (fl"3j) is also shown. We observe that the prediction 
•AeffjC = 1 ° n ly holds for low values A^Kem while for 
larger values the simple requirement N c ff > 1 is not 
enough. A careful study revealed that for large values 
of Ti s i n and iii n the critical effective soliton number 
scales as Af e g iC oc Ti^ii^. This suggests that the ef- 
fect is not related to any of the soliton numbers (which 
scale as Ti^ ny yii/ m ). Instead, consider the energy flu- 
ence <f>j = J dtlj(t), which for a sech-shaped pulse 
Iosech(i/T ) is $ = 2T I . In Fig. [3] the corresponding 
N s, c from the Afc c data of Fig. [5] is plotted vs. 
and the data are well represented by the rational fit 



AU c = 1 



where the free fitting parameter 4> 



*' 1)in = 4>i,i„/$ c , (15) 
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Fig. 3. The data from Fig. [2] shown vs. the input energy 
fluence $i,i n = 2Xij n /i, in in a log-log plot. 



expect that this <I> C value is only valid for the partic- 
ular case studied. A possible explanation for the devi- 
ation from N e s lC = 1 was found by studying the ini- 
tial stages of propagation when the fluence is high such 
that 1 < AT cff < N oStC , see Fig. H for Afc = 137 mm" 1 
and $1,111 = 272 mJ/cm 2 . As one would expect when 
N c g > 1 a substantial negative FW phase shift builds 
up. However, in the pulse center an increase in the phase 
is seen, which causes the chirp across the pulse center 
to be non-monotonic making the pulse unable to com- 
press upon further propagation. By turning off the Kerr 
cross-phase modulation (XPM) terms in the numerics, 
we observe that this phase increase in the pulse center 
disappears; now the pulse has a monotonic chirp across 
the pulse center so it can compress. This was indeed 
observed upon further propagation. Turning on the Kerr 
XPM terms again and taking Afc = 135 mm -1 to achieve 
N e ff > N D ff c ^> 1, N c ff is now large enough cancel also 
the detrimental contribution from the Kerr XPM terms: 
the phase increase in the center is no longer there mak- 
ing the chirp monotonic. Therefore the pulse can com- 
press, which we indeed observed. Thus, the deviation 
from N e g lC = 1 comes from a positive phase contribu- 
tion in the pulse center, which seems to originate from 
the Kerr XPM terms. These cannot be neglected due to 
the high fluence, but were indeed neglected in the NLSE- 
type model ([9]) from which the conjecture N e s tC — 1 orig- 
inated. An indication of this high-fluence effect on the 
XPM terms was recently observed experimentally [23]. 
Finally, as the Kerr XPM coefficient is changed we found 
that $ c oc B- 1 / 2 . 

We can conclude that N c g, c deviates quite strongly 
from unity when the fluence (and thus the soliton num- 
bers) is large. However, for $i,i n < 3>e, Aeff.c = 1 can be 
used MS IX safe estimate of the transition point . 

B. Scaling laws for compression parameters 

We proceed to investigate how the soliton compressor 
behaves when compression is successful. We seek to find 
general scaling laws that can tell us at what position 
z op t in the crystal the pulse reaches its optimal single- 
spike compressed state - i.e., with maximum intensity 
and minimum duration - and to check its duration, com- 
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Fig. 4. Numerical simulations showing the FW phase and 
chirp after propagation in 50 mm BBO with a high flu- 
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pression factor, and quality. 

For soliton compression in the NLSE Mollenauer et 
al. [12] first studied the compression for 1 < A^err < 15. 
Based on this Dianov et al. [14] showed the following 
empirical scaling law for the optimal compression length 



^opt 

zo 



0.32 



Kerr 



l.i 



(16) 



where z$ = §£d.i is the soliton length [10]. Also based 
on Ref. [12], Tomlinson et al. [13] reported the empirical 
scaling law for the compression factor f c = Ti t i n /Ti^ opt 



f c = 4.1N KeiI , 1<AW<50. 



(17) 



The validity range was reported in Ref. [14]. Finally, 
the pulse quality is defined as the fluence in a sech 2 - 
fit to the central spike relative to the input fluence 
Q c = $i, 8 ech— fit/^i,in> which is equivalent to taking 
Q c = ii/Ji^n/c- Thus, we should expect that Q c oc l// c 
[11]. As we will now show, the cascaded quadratic soli- 
ton compressor follows nicely these scaling laws when 
A^Kerr is replaced by AT c ff. This holds as long as we are 
in the stationary regime, i.e., when the NLSE-likc model 
Eq. ((9j) is a good approximation to the system and when 
GVM effects are small. 

We performed a wide range of simulations using sech- 
shaped un-chirped input pulses with parameters ranging 
from T Un = 50 - 2000 fs and 7 Un = 10 - 400 GW/cm 2 . 
The simulations used in the following resulted in cleanly 
compressed pulses as demonstrated in Fig. [5] It is worth 
to stress that the depletion of the FW due to conversion 
into SH was quite low in all cases (< 3% in the station- 
ary regime) due to the large values of phase mismatch 
40 < Afc < 120 accessible in the compression window. 
Some selected compression examples are given in Fig. [5) 
(a) shows a typical case with a 200 fs input, where we 
have optimized the phase mismatch and the input in- 
tensity to give a 6.0 fs compressed pulse. This implies 
a compression ratio of f c = 33 and the pulse quality is 
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Fig. 5. Selected simulations of clean compressed pulses 
at the optimal compression length, (a) Afc = 50 mm , 
Ti )in = 200 fs FWHM and 7i, in = 59 GW/cm 2 (N eS = 8) 
resulting in a 6.0 fs pulse (/ c = 33) with Q c = 0.26. (b) 
Afc = 55 mm- 1 , T Un = 2000 fs FWHM and I 1M = 
26.7 GW/cm 2 (A^ eff '= 50) resulting in a 7.6 fs pulse 
(f c = 264) with Q c = 0.06. 

Q c = 0.26. The remaining pulse energy resides in the 
unwanted pedestal as well as in the SH (2.3% conver- 
sion occurred). As a more extreme case, (b) shows a 2 
ps long pulse compressed to a clean 8.3 fs pulse, imply- 
ing an impressive f c = 264. A lot of the energy remains 
in the pedestal (only 0.5% is converted to the SH), and 
the pulse quality is merely Q c = 0.06. This is typical 
of large compression ratios. Note the difference in z opt ; 
z pt = 342 mm in Fig. [5](b) is long, but may be realized 
using multiple crystals. 

In Fig.^a) we show z opt /z as function of the effective 
soliton number N e s . The data follow the scaling law (p~6|) 
- with A^Korr is replaced by N e e - quite well, but there 
are deviations for small N c ff. This is due to the A^~ 2 
term in Eq. (|16|). so an improved fit gave 



^opt 



0.44 2.56 



0.002. 



(18) 



We have also included data points where the compressed 
pulses were less clean (cither having trailing or leading 
oscillations, or being somewhat asymmetric). As the plot 
indicates they start kicking in when A^g > 10, which can 
be explained by increased XPM contributions as well 
as increased influence of the GVM induced Raman-like 
perturbation. Nonetheless, these less clean data still fol- 
low the scaling law (fT8|) closely. An exception is when 
the nonstationary regime is close: We included four data 
points (black symbols), which have the same parame- 
ters as Fig. EJa) and where Afc is gradually decreased 
so the final data point is in the nonstationary regime. 
These data points start to deviate from the scaling laws 
because the compressed pulses experience more strongly 
the GVM induced Raman-like effects. 

In Fig. [Hb) we plot the compression factor f c as func- 
tion of AT c ff. Again the scaling law (fT7|) holds well except 
for small N e g. A linear fit to the clean data gave the 
following scaling law valid also for small A^ e ff 



/ c = 4.7(AT cff -0.86). 



(19) 
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Fig. 6. Results of numerical simulations showing the op- 
timum compression parameters vs. N e f{ in log-log plots 
for (a) the compression length z op t/zo, (b) the compres- 
sion factor f c and (c) the compression quality Q c . The 
simulations marked with round symbols resulted in a 
clean compressed pulse, while the triangles resulted in 
less clean pulses. The four black triangles gradually en- 
ter the nonstationary regime. The solid lines are fits to 
the clean data [Eqs. (p~5|) - (|^T7j) ] . while the dotted lines are 
the scaling laws (|T6|) and (fl7|) . 



The less clean data also follow the scaling law (flU)) quite 
well, but the data points approaching the nonstationary 
regime (black triangles) very quickly separate out. Thus, 
f c is very sensitive to the effects of the GVM induced 
Raman-like effects in the nonstationary regime. 

As the last dimensionless parameter, we show in 
Fig. [6jc) the pulse compression quality Q c as a function 
of N e ff . The data roughly follow a rational function 



Q c = [0.24(7V cff - l) 1 



(20) 



The exponent 1.11 deviates from unity due to the be- 
haviour of Q c for small N e s. However, for large N e g we 
find Qc oc /~ , as predicted. 

As a more concrete example Fig. [7] shows the com- 
pressor length one should choose for optimal compres- 
sion z opt , together with the expected compressed pulse 
duration. These were calculated from Eqs. (fT8"|) and (fT9")l 
for Afc = 50 mm -1 . Single-cycle pulses are available for 
T^in > 200 GW/cm 2 with realistic compressor lengths 
around 5 — 20 mm. 

Finally, Fig. [8] summarizes output pulse duration vs. 
N e s, providing a sense of the possible pulse compression 




Fig. 7. The optimal compressor length and the expected 
compressed pulse duration in a BBO for Ai = 1064 nm 
and Afc = 50 mm -1 with (a) Ti^n fixed; (b) Ix,xn fixed. 
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Fig. 8. The duration of the compressed pulse vs. N e s in 
a log plot. The data are from the same simulations as in 
Fig.[6l but only those with Afc = 55 mm -1 are shown. 



in BBO at Ai = 1064 nm. The lowest value observed 
was 4.7 fs FWHM, corresponding to 1.3 optical cycles. 
Single-cycle pulses were not observed, most likely due to 
GVM-induced Raman effects and Kerr XPM effects. 

4. Conclusions 



In summary, we find that the effective soliton num- 
ber JVcff = y/ A^ht: — -^Kcrr ^ s the P r0 P cr dimensionless 



IV SHG 

quantity for describing soliton compression using cas- 
caded quadratic nonlinearities in the stationary regime. 
Soliton compression generally only occurs when 7V ff > 1, 
and this in turn requires that the quadratic nonlinearity 
dominates over the cubic nonlinearity. In this balance, 
material and input pulse parameters such as the phase 
mismatch, the input pulse intensity and duration play 
a crucial role. When the pulse energy is large N e g > 1 
is no longer a sufficient demand. We attributed this to 
Kerr XPM effects, and we found an empirical scaling law 
relating N e s tC to the input fluence. 

We showed through a large number of realistic numeri- 
cal simulations of a BBO crystal that the system clearly 
obeys dimensionless scaling laws dictating the optimal 
compression propagation distance, the pulse compres- 
sion factor and quality. These scaling laws were expressed 
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through the effective soliton number iV G ff, and are very 
similar to the ones observed in the NLSE. The scaling 
laws are general and hold also for other materials and 
wavelengths, and will serve as a crucial tool to deter- 
mine input pulse peak intensity and duration as well as 
phase-matching conditions and optimal crystal lengths. 

Besides the requirement of a strong quadratic nonlin- 
carity, avoiding the nonstationary regime, where GVM- 
effects are strong, is another obstacle to observe the de- 
sired compression in the considered system. While com- 
pression may still occur, the final pulse is too distorted 
to be of any use due to GVM induced Raman-like effects. 
This is particularly a problem for short pulses, high in- 
tensity pulses, and in materials with a weak quadratic 
nonlincarity and/or large GVM. In this paper we consid- 
ered the stationary regime, where this effect is negligible 
and stress that the scaling laws for compression hold only 
in this regime. In subsequent publications we will study 
the GVM effects in more detail, and show that the cas- 
caded quadratic nonlincarity induces a nonlocal response 
on the FW, and that in presence of GVM this response 
becomes Raman-like [16]. Moreover, in the nonstation- 
ary regime the response function is no longer localized, 
which has severe implications on the compressed pulse. 

Let us finally touch on the exciting prospect of the 
cascaded quadratic soliton compressor in a wave-guiding 
configuration. As mentioned, GVM effects are detrimen- 
tal to the compression performance, and these become 
even more pronounced for shorter input wavelengths 
such as 800 nm because GVM is much stronger. How- 
ever, in a photonic crystal fiber (PCF) wave-guiding ef- 
fects can dramatically alter the dispersion, and in partic- 
ular it was shown that for SHG it is possible to achieve 
zero GVM [24]. Moreover, in a PCF the mode areas are 
very small (roughly a few /im 2 is possible) so very high 
peak intensities are possible even with very low pulse en- 
ergies. This will allow achievement of low-energy clean 
compressed pulses, once issues with quadratic nonlinear 
response of optical fibers have been resolved. 
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A. The cubic nonlinear response 

We here show how a cubic nonlinearity is included in the 
model. The cubic nonlinear polarization response is 



p(3) _ . 



dtidt 2 dt^ 



x x} >(t -t u t- t 2 ,t - i3):E(*i)E(ia)E(t 3 ), (21) 

%( 3 ) is a rank 4 tensor describing the cubic nonlinear 
response of the material, which we let take the func- 
tional form x (3) (ii^2,t3) = X^ R(h)S(t 2 )8(t 3 ), where 



R(t) is the normalized material Kerr response func- 
tion. In arriving at Eqs. ([T]) we then assumed the waves 
monochromatic and polarized along arbitrary directions 
E(<) = Re^-Ea^e-^ 1 * + x 2 E 2 (t)e~ iU2t ], where Xj is 
the unit polarization vector. In SHG the FW and SH 
waves are either all polarized along the same direction 
(Type 0), or the two FW are polarized along the same di- 
rection but orthogonal to the SH (Type I), or finally the 
two FW photons could be polarized orthogonally to each 
other while the SH is parallel with one of the FW pho- 
tons (Type II) [19]. The propagation equations presented 
here are valid for Type and I, since the FW photons 
are considered degenerate. We assume an isotropic cubic 
Kerr nonlinearity, for which the only nonzero tensor com- 
ponents of x (3) are Xxxyy = Xxyyx = Xxyxy = Xxxxx/3 
[11]. This means that the XPM coefficient in Eq. ([Ta|) 
is B = 2 for Type SHG, and B = 2/3 for Type I 
SHG. We now divide the material Kerr response into an 
electronic response and a vibrational Raman response 
R(t) - (1 - f R )S(t) + f R h a (t) [10], where f R is the frac- 
tional contribution of the vibrational Raman response. 
The first part describes the electronic response, which 
can be considered instantaneous, resulting in the cubic 
SPM and XPM terms in Eqs. (HJ). The vibrational Ra- 
man response for a general 3-wave mixing process with 
two different frequencies is [25] 



Hj{r) = Sj / dsh R (s) 

J — oo 

x{E ] (r)[\E J (r- S )\ 2 + ^B\E m (r- S )\ 2 ] 
+ \BE 3 {r - s)E* m (r - s )e^~^ s E m (r)} , (22) 

where j = 1,2 and m = 3 — j, and only terms that are 
phase matched are included [25]. This Raman response 
is governed by the gain function h R . For SHG, the in- 
tegral over the term containing e *("j- w ™) s i s vanishing 
since it implies evaluation of h R is evaluated at frequency 
offset \ujj — uj m \ = uj\ much larger than the typical spec- 
tral width of h R (in the THz range). Therefore, for SHG 
only the term Ej(r) [\Ej (r - s)| 2 + \B\E m (r - s)\ 2 ] re- 
mains, and the generalized coupled equations (|TJ) become 



C\Ei + Kg HG1 SiElE 2 e lAkz 
+ 4 er r,i&£i | (1 - /*) (\Ei\ 2 + B\E 2 \ 2 ) 



(23a) 



f R ^ dsh R (s) [{E^t -s)\ 2 + \B\E 2 {t - s)| 2 ] 1 

J —oo J 



2 fi 2 + ^SHG. 2*2-^16 



= o, 

(23b) 



K^S.EJil-f^HE^+BlE^) 



f R J dsh R (s) [|£i(t -s)\ 2 + \B\E 2 {t - s)| 2 ] | = 0. 
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Finally, the Raman overlap integral has an approxi- 
mate form for short, but not extremely short, pulses 

f B K 3 {T) ~ f R S 3 E 3 (r) [\E 3 {r)\ 2 + \B\E m {r)\ 2 ] 

-r R -^[\E 3 (r)\ 2 + ^B\E m (r)\ 2 ], (24) 

where we have used that J^dth^t) = 1, and de- 
fined the well-known Raman response time t r = 
fn I^oa dtth R (t). Eq. is related to intra- and inter- 
pulse Raman scattering (IIRS), and Eqs. ([23]) become 

CiEk + Ki liGA S 1 E*E 2 e tAkz (25a) 

+ 4err,l{£li?l [|i?l| 2 + B(l - \f R )\E 2 \ 2 ] 

-r il S 1 ^ : (| J B 1 | 2 + i J B|£; 2 | 2 )} = ) 
C 2 E 2 + K i Wi2 S 2 E 2 ie - iAkz (25b) 
+ «K™r, a {^^ [\E 2 \ 2 + B(l - y R )\E x \ 2 ] 

- trE^^E^ 2 + \B\E 2 \ 2 )} =0. 

Here we have neglected the 2nd order derivatives from 
applying self-steepening to the IIRS term in Eq. (|2"4"f . 

B. Extending the propagation equations to the 
SEWA regime 

Here we discuss the difference between the slowly varying 
envelope approximation (SVEA) and the SEWA propa- 
gation equations. SEWA [17] is a general spatio-temporal 
model that describes spatio-temporal pulse propagation 
down to single-cycle pulse durations. It was recently de- 
rived for SHG by Moses and Wise [18]. SEWA does not 
pose any direct constriction on the pulse bandwidth, 
whereas the extended SVEA (with steepening terms and 
the general Raman convolution response) holds for [26] 
Aoj/uj < 1/3. Neglecting transverse spatial terms, the 
only term the SHG extended SVEA model does not in- 
clude is related to the dispersion of the SH, since due 
to GVM the SH dispersion operator ^ must be re- 
placed by the following effective operator [18] T> 2 , fi = 

i* 2 + S^ 1 ^^-^. Imposing the scalings z' = z/L-qi and 
t' = T/Ti.i n , the dimcnsionlcss operator becomes 

% cS ^V' 2 + S' 2 ' V -^, (26) 

where the dimensionless factor v = cd 2 2 /uj 2 n 2 \k^\. 
Expanding the inverse steepening operator S!^ 1 = 

=£.) ^r, we get 

%* = £ >'•" 4"" + 5 (|) («) 

m=2 L V 7 J 

In the SEWA there is no restriction on the pulse band- 
width so single-cycle temporal resolution is achieved. 



However, since we consider SHG we should be careful. 
One assumption made when deriving Eqs. |T]) is namely 
that the spectra of the FW and SH do not overlap (sub- 
stantially). This assumption allows us to separate the 
fields as shown in App. [A"l We chose Aui/uij = 0.7. This 
gives some overlap between the FW and SH spectra, but 
we always made sure that the spectral components in 
the overlapping regions were negligible. 

C. BBO 

We use a (3-BaB 2 C>4 (beta-barium-borate, BBO) non- 
linear crystal, where collinear Type I SHG is possible 
through the interaction oo — > e, i.e. the FW photons 
are ordinarily polarized, while the generated SH photon 
is extraordinarily polarized. BBO is a uni-axial crystal 
where phase matching can be achieved by bircfringent 
phase matching by changing the angle 9 between the 
FW input and the optical z-axis of the crystal. d e ff also 
changes with 8, and is on average 2.22 pm/V in the area 
in which we are interested in. The dispersion is calculated 
from the Sellmeier equations [19]. Note in this connec- 
tion that the analytical transition to compression V e ff jC , 
as calculated from Eq. (|13|) . is on implicit form. 

BBO is an excellent nonlinear medium for the cur- 
rent purpose because it has a decent quadratic nonlinear 
strength, and perhaps more important, it has a small 
GVM at NIR wavelengths, which is the major reason for 
using BBO instead of, e.g., periodically poled LiNb03 
for these wavelengths. BBO also has a very low cubic 
nonlinear refractive index. As we noted from Eq. (|10p 
the balance between these is crucial. BBO also has a 
very low two-photon absorption coefficient (except in the 
ultra-violet part of the spectrum [27]), justifying the ap- 
proximation made in the derivation of Eqs. ([S"|). 

In the literature several values for the nonlinear re- 
fractive index have been reported. We chose to use 
n^ e „ A = 3.65 ± 0.6 ■ 10~ 20 m 2 /W reported in Rcf. [28] 
mainly because the measurements are done at phase- 
matching angles close to ours and also used fs pulses, 
but also because this value in the past has given the 
best agreement with experimental observations. 
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